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Abstract 

We investigate the conditions for reactivity enhancement of catalytic processes in porous solids 
by use of molecular traffic control (MTC). With dynamic Monte-Carlo simulations and continuous- 
time master equation theory applied to the high concentration regime we obtain a quantitative 
description of the MTC effect for a network of intersecting single-file channels in a wide range 
of grain parameters and for optimal external operating conditions. Implementing the concept of 
MTC in models with specially designed alternating bimodal channels we find the efficiency ratio 
(compared with a topologically and structurally similar reference system without MTC) to be 
enhanced with increasing grain diameter, a property verified for the first time for an MTC system. 
Even for short intersection channels, MTC leads to a reactivity enhancement of up to approximately 
65%. This suggests that MTC may significantly enhance the efficiency of a catalytic process for 
small as well as large porous particles with a suitably chosen binary channel topology. 
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I. INTRODUCTION 



Zeolites are used for catalytic processes in a variety of applications, e.g. cracking of large 
hydrocarbon molecules. In a number of zeolites diffusive transport occurs alongquasi-one- 
dimensional channels which do not allow guest molecules to pass each other Due to 
mutual blockage of reactants A and product molecules B under such single-file conditions 
P| the effective reactivity of a catalytic process A — ► B - determined by the residence time 
of molecules in the zeolite - may be considerably reduced as compared to the reactivity 
in the absence of single-file behaviour. It has been suggested that the single-file effect 
may be circumvented by the concept of molecular traffic control (MTC) 0, 0] which has 
remained controversial for a long time. This notion rests on the assumption that reactants 
and product molecules respectively may prefer spatially separated diffusion pathways and 
thus avoid mutual suppression of self-diffusion inside the grain channels. 

The necessary (but not sufficient) requirement for the MTC effect, a channel selectivity 
of two different species of molecules, has been verified by means of molecular dynamic (MD) 
simulations of two-component mixtures in the zeolite ZSM-5 [f| and relaxation simulations 
of a mixture of differently sized molecules (Xe and SF 6 ) in a bimodal structure possessing 
dual-sized pores (Boggsite with 10-ring and 12-ring pores) (a). Also equilibrium Monte-Carlo 
simulations demonstrate that the residence probability in different areas of the intracrys- 
talline pore space may be notably different for the two components of a binary mixture [3] 
and thus provide further support for the notion of channel selectivity in suitable bimodal 
channel topologies. 

Whether an MTC effect leading to reactivity enhancement actually takes place was ad- 
dressed in a series of dynamic Monte Carlo simulations (DMCS) of a stochastic model system 
with a network of perpendicular sets of bimodal intersectingchannels and with catalytic sites 
located at the intersecting pores (NBK model Fig. [T] b)) Q, S H Q Q. The authors of 
these studies found numerically the occurrence of the MTC effect by comparing the outflow 
of reaction products in the MTC system with the outflow from a reference system (REF 
model Fig. da)) with equal internal and external system parameters, as well as equal internal 
and optimal external system parameters. The dependency of the MTC effect as a function 
of the system size and catalytic reactivity has been investigated in (ll|, The MTC effect 
is favored by a small number of channels and long channels between intersections, which by 
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FIG. 1: REF system a), NBK model consisting of perpendicular a//3-channels b) - both with 
channel number of N = 3. In contrast to the REF case, where we allow both types of particles 
(A and B particles) to enter any channel, in the MTC system A particles are carried through the 
a- channels whereas the B particles diffuse along the (3- channels. Black squares indicate catalytic 
sites where a catalytic transition A — » B takes place. 

themselves lead to a very low absolute outflow compared to a similar system with shorter 
channels. For reasonable reactivities and channel lengths the MTC effect vanishes inversely 
proportional to the grain diameter. An analytical treatment of the master equation for this 
stochastic many-particle model revealed the origin of this effect at high reactivities [l^j]. It 
results from an interplay of the long residence time of guest molecules under single-file con- 
ditions with a saturation effect that leads to a depletion of the bulk of the crystallite. The 
extension of the NBK model to three dimensions |14] leads to similar results with the con- 
clusion that MTC may enhance significantly the effective reactivity in zeolitic nanoparticles 
with suitable binary channel systems and thus may be of practical relevance in applications. 
However, no MTC effect with sufficient absolute effective reactivity can be expected with 
NBK channel topology for large grains of at present commercially used size. 

Here we address this question by investigating the MTC effect for a novel channel-topology 
(Fig. EJ) consisting of alternating a-and /3-channels. With a view on applications we focus 
on small channel length. Moreover, for the same reason we determine the MTC effect 
by making a comparison with the reference system using the same set of fixed internal 
(material-dependent) parameters, but (unlike in previous studies [ill [ll|, [isj]) for each case 
(MTC and REF resp.) different optimal external (operation-dependent) parameters which 
one would try to implement in an industrially relevant process. This may be of interest as 
since the first successful synthesis of mesoporous MCM-41 nanoparticles there has been 




FIG. 2: MTC system with alternating a//3-channels (N = 5). 

intense research activity in the design and synthesis of structured mesoporous solids with a 
controlled pore size. In particular, synthesis of bimodal nanostructures with independently 
controlled small and large mesopore sizes has become feasible It will transpire that a 
signifcant MTC effect (reactivity enhancement up to 65%) occurs in our model system even 
for small channel length at realistic intermediate reaction rates of the catalyst. 

II. THE MODEL 

We consider a lattice model with a quadratic array ofNxN channels (Fig. [2} which is a 
measure of the grain size of the crystallite. Each channel has L sites between the intersection 
points where the irreversible catalytic process A ^ B takes place. We assume the boundary 
channels of the grain to be connected to the surrounding gas phase, modelled by extra 
reservoir sites of constant A particle density p. We assume the reaction products B which 
leave the crystallite to be removed immediately from the gas phase such that the density 
of product particles in the reservoir is always 0. Short-range interaction between particles 
inside the narrow pores is described by an idealized hard-core repulsion which forbids double 
occupancy of lattice sites. 

The underlying dynamics are stochastic. We work in a continuous-time description where 
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the transition probabilities become transition rates and no multiple transitions occur at the 
same infinitesimal time unit. Each elementary transition between microscopic configurations 
of the system takes place randomly with an exponential waiting-time distribution. Diffu- 
sion is modelled by jump processes between neighbouring lattice sites. D is the elementary 
(attempt) rate of hopping and is assumed to be the same for both species A, B of particles. 
In the absence of other particles D is the self-diffusion coefficient for the mean-square dis- 
placement along a channel. If a neighboring site is occupied by a particle then a hopping 
attempt is rejected (single-file effect). The dynamics inside a channel are thus given by the 
symmetric exclusion process which is well-studied in the probabilistic Q 

and statistical mechanics literature |22|. The self-diffusion along a channel is anomalous, 
the effective diffusion rate between intersection points decays asymptotically as 1/L, see [3| 
and references therein. 

At the intersections the reaction A — > B occurs with a reaction rate c. This reaction rate 
influences, but is distinct from, the effective grain reactivity which is largely determined 
by the residence time of guest molecules inside the grain which under single-file conditions 
grows in the reference system with the third power of the channel length L j^J. At the 
boundary sites, particles jump into the reservoir with a rate D(l — Pa~ Pb) in the general 
case. Correspondingly, particles are injected into the grain with rates Dpa,b respectively. 
As discussed above here we consider only pa = p, Pb = 0. 

For the REF system reactants as well as product particles are allowed to enter and 
leave both types of channels. In case of MTC reactant (product) particles will enter a(/3)- 
channels only, mimicking complete channel selectivity. Therefore all channel segments carry 
only one type of particles in the MTC case. For the boundary channels complete selectivity 
implies that a-channels are effectively described by connection with an A-reservoir of density 
Pa = p (particles of type B do not block the boundary sites of a-channels) and /3-channels 
are effectively described by connection with a S-reservoir of density ps = 0, respectively, 
(particles of type A do not block the boundary sites of /3-channels.) This stochastic dynamics, 
which is a Markov process, fully defines our MTC model. 

In both cases, MTC and REF, the external concentration gradient between reactant and 
product reservoir densities induces a particle current inside the grain which drives the system 
in a stationary nonequilibrium state. For this reason there is no Gibbs measure and equi- 
librium Monte-Carlo algorithms cannot be applied for determining steady-state properties. 
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a) b) 

Profile of A particles - MTC2 (N=L=9, c=0.1) Profile of B particles - MTC2 (N=L=9, c=0.1) 




FIG. 3: A/B particle profile (a)/b)) for the MTC system of N = L = 9 and reactivity c = 0.1. 

Instead we use dynamic Monte-Carlo simulation (DMCS) with random sequential update. 
This ensures that the simulation algorithm yields the correct stationary distribution of the 
model. 



III. MONTE CARLO RESULTS 



Anticipating concentration gradients between intersection points we expect due to the 
exclusion dynamics linear density profiles within the channel segments [isl [li], [3] ? the slope 
and hence the current being inversely proportional to the number of lattice sites L. The 
total output current j of the product particles, defined as the number of B particles leaving 
the grain per time unit in the stationary state, is the main quantity of interest. It determines 
the effective reactivity of the grain. 

Fig. El shows the particle densities for the reactants a) and product particles b) in the 
stationary regime for the MTC system N = L = 9 and reactivity c = 0.1. Every other 
boundary-channel segment contributes to the output, therefore we expect no efficiency loss 
with increasing N. We are particularly interested in studying the system in its maximal 
current state for given reactivity c and size constants N, L, which are intrinsic material 
properties of the grain. The reactant particle reservoir density p, determined by the density 
in the gas phase, can be tuned in a possible experimental setting. Let us therefore denote 
the reservoir density which maximizes the output current with p* and the maximal current 
with j*. We iteratively approach the maximal current by a set of 9 datapoints. Fig. |U 
shows j as a function of p for both a MTC and REF system of N = 5, L = 1 and reactivity 
c = 0.01. The "best" datapoint has been chosen in order to approximate the maximum. 
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FIG. 4: Jref (solid symbols) and Jmtc (open symbols) as a function of the reservoir density for a 
system with N = 5, L = 1, c = 0.1. 
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j for L=2 and N=3.. 15 (total selectivity) 













N=15 




2 
1.5 

> 1 












' N=1 5 




















































> 


• * 












e 

e 









n 








* 


> » 






I, 


I \ I 


e 

B B 


B ° 


3 B B 


N=3 




. i, 


1 i.S 


B H 


B B 


B B 


N=3 





0.0001 0.001 0.01 0.1 1 10 100 1000 0.0001 0.001 0.01 0.1 1 10 100 1000 

c c 

FIG. 5: Maximal currents j* over reactivity c for the REF system a) and MTC system b) for 
different N. 

Statistical errors are displayed, they are, however, mostly within symbol size. We note that 
for the MTC system the maximal current is attuned for p* — 1. REF and MTC system 
show a similar dependence of the maximal current as a function of the reactivity (Fig. [HJ). 
In contrast to the NBK topology [13], j*(c) is a monotonically increasing function and the 
highest product particle output is achieved for infinite reactivity. In this case any transition 
A^B transpires close to the surface and hence, from an applied point of view one would be 
interested in various small grains rather than view big ones. 

In order to measure the efficiency of an MTC system over the associated REF system we 
define the efficiency ratio 



R(c,N,L) = 3 -MIc 
Jref 



(1) 



which is a function of the system size N, L and reactivity c. 

Fig. El a) shows the measured ratio R for a large range of reactivities. We plotted systems 



j MTC /j REF for L=2 and N=3..15 (total selectivity) 



Maximal Ratio j MTC /j REF L=1,2,3 (total selectivity) 
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FIG. 6: a) Ratio R(c) for different number of channels N and L = 2. b) Maximal ratio R* for 
different L 

with L = 2 and different N. We find an MTC effect over a large range of reactivities. There 
is a small tendency of growing efficiency with increasing N. This is highly desirable since 
an enhancement of efficiency due to MTC is thus not restricted to nano-scale crystallites. 
Compared to former studies of MTC systems, this is a qualitatively new property and results 
from the fact that the number of output (/3-channels) scales linearly with growing system 
size N. In general, compared to the NBK model, the efficiency ratio keeps quite unchanged 
in the regime of simulated reactivities. There is no sharp decrease for large reactivities. 
Nevertheless, we notice an optimal reactivity c* for which R becomes maximal. This value 
is denoted by R* and plotted in Fig. Elb) for increasing N and different L. The MTC effect 
is present for any L and our model exhibits an enhancement of the effective reactivity of up 
to 20% for (extremely) short channel lengt h ( L = 1) and up to 65% for L = 3. 

As expected from previous results |ll|, ll2j the MTC effect is seen to increase with in- 
creasing L (Fig. EJ) . This follows from theoretical studies of single-file systems which predict 
that the mean travelling time of a product molecule through a channel of length L is pro- 
portional to L 2 in the MTC case as in ordinary diffusion, but proportional to L 3 in the REF 
case due to mutual blockage. Hence the current is proportional to 1/L in an MTC system, 
but proportional to 1/L 2 in a REF system. This holds for all values of the parameters and 
hence, for sufficiently large L, the MTC system becomes more efficient. 
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j/j REF for L=1 ,2,3 and N=15 (total selectivity) 
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FIG. 7: Dependence of the output ratio on the channel length L. 




FIG. 8: A/B particle profile (a)/b)) for a system of N = 9, L = 15 and reactivity c = 1000. 
IV. THE HIGH-REACTIVITY CASE 

The penetration depth of reactant particles into the system is related to the reactivity 
c. In the fast-reactivity case transitions and particle transport take place close to the grain- 
surface. Fig. [HI shows the profile of the reactant particles a) and product particles b) for a 
system of iV = 9 channels, L = 15 and fast reactivity (c = 1000). Notice that almost no 
reactants reach the bulk of the system since reactants entering the system are converted very 
likely into product particles at the first intersection. Hence, a creation of product particles 
occurs only on the outer /^-channels, which leads to equilibration of the product particles 
in the bulk of the system. This picture is supported by simulation (Fig. (HJ) and suggests to 
define a reduced system, illustrated in Fig. EE for the following theoretical analysis. 

Here, a-segments carrying no particles are suppressed, /9-segments in equilibrium are 
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FIG. 9: Corner of the reduced system for the case of fast reactivity. 

represented by dotted lines. The density profile in a channel segment enclosed by two 
intersections or an intersection and boundary site is linear due to pure exclusion dynamics 
This allows us to reduce our effort to the computation of the intersection densities. 
The intersection densities are calculate from the Master equation. The output of product 
particle, i.e. the number of B's .saving the system, is given by the current, g of the 
/5-segments connected to the reservoir. The calculation is similar to [12] and we identify 

Here, p x is the product particle density at intersection x connected to the /^-Segment. The 
derivation of the currents involves a simplification. Within a mean field approximation, 
joint probabilities (xy) of particle numbers between an intersection and its adjacent sites 
can be replaced by the product of the individual averages (x) (y). This well controlled 
approximation is, in fact, necessary only in the vicinity of intersections, since any joint 
probability vanishes in the bulk of the channels 1 201 . The expected relative error is of the 

nn 

order of 1/L |U2, |22j . Reactant particles diffuse into the system through a-segments which 
are bounded at one side by the reservoir (density pa) and at the other side by intersection 
p x with current 

a _ D Pa(Px-I) (3) 

3x D L{p x -i)-r {3} 
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Finally, the intrinsic currents, determined by two adjacent intersections follow to 

7S = D P ~j^ (4) 

L + l w 

The arrows are not supposed to indicate vectors but to distinguishe the direction of currents 
in order to comply with the zigzag profile shown in FigJHl 
Making use of conservation of currents 

m = 7: + 7: (6) 

Jx 7" ■ 7" :• (7) 

we find (for large but finite L) two coupled recursion relations for the intersection densities 
which determines the density profile. 

Jl 

P* = \ (8) 
I \(p x + Px+i + Pa) 

3px = px + Px-i (9) 

Let N a = \_N/2\ be the number of a-channels. In order to ensure zero current in the 
equilibrium parts of the /3-channels (Fig. EJ, intersections marked with a black square need 
to attain equal densitiy. Hence, in the region 2 < x < N a we assume a constant solution 
for the densities p x and find p x = p = pa- Intersections p x connected to a-channels carry 
a product particle density being 3/2 times higher (Fig. HHJ and therefore p x = p = 3/2pA- 
Since the particle density can not exceed 1, p saturates for a boundary density of pA = 2/3. 
This has an interesting consequence. Changing the boundary density from pA = 2/3 to 
Pa — 1 merely shifts up the profile of the a-channels, the output, however, keeps unchanged. 
For the sake of completeness we remark that for the four corners 2p\ = p 1? hence, 3pi = 2p^. 

We have now determined the profile. Because of (J2J) the output of product particles for 
each /^-channel is independent of the number of channels, hence, the total output current 
scales linearly with N. Figure HH a) shows the profile of intersection p 2 as a function of 
the reservoir density for a system of iV = 9 channels and different L. The saturation at 
Pa = 2/3 predicted for L — > oo (vanishing corrections to mean field theory) becomes more 
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FIG. 10: Density p2 for different L as a function of the reservoir density. The theoretical density is 
indicated by the dotted line. N = 9 

pronounced for increasing L. The theoretical description is good for small reservoir densities. 
The largest deviation (about 20% at most) occurs for the predicted saturation point, so for 
practical purposes the theoretical analysis agrees well with the simulated data. 



V. MTC MODEL WITH SIZE SELECTIVITY 



For situations where channel selectivity merely originates from different pore diameters 
we have to modify our model. Consider the model shown in Fig. [T] but now, allowing 
reactant as well as product particles to diffuse along the a-channels (size selectivity). We 
call them 7-channels. Assuming /3-channels to be smaller in diameter than a-channels there 
is no reason why particles of type B would avoid to occupy a-channels. Of course we could 
also consider the situation vice versa, with (3 being larger than a. However, this leads to 
a reduction of output channels, hence, this situation would be counterproductive and is 
not studied here. Overtaking is still prohibited for any type of particles. The dynamics 
explained in the model-section keeps unchanged and similar to arguments given before we 
would like to compare with the REF system for equal internal parameters (N, L, c) but in 
their maximal current state with respect to the reservoir density pa- 

Comparing the currents of the size-controlled system (Fig. [TT] a)) with the model of 
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j for L=2 and N=3..15 (size selectivity) 



j for L=1 ,2,3 and N=1 5 (total/size selectivity) 
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FIG. 11: a) Maximal currents j* over reactivity c for the size controlled MTC system, b) Maximal 
Currents for the size controlled system (solid symbols) and the system of total selectivity (open 
symbols). 
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FIG. 12: Left: Ratio -R(c) for different number of channels N and L = 2. b) Maximal ratio R* for 
different L. (size selectivity) 

total selectivity (Fig. El a)) we find a similar qualitative behaviour. There is an interesting 
observation, that the size-controlled model performes a higher output compared to the model 
of total selectivity for L = 1, however, for larger L total selectivity becomes more and more 
efficient (Fig. [H]b)). This is true for high enough reactivities. Both models maximize their 
output for infinite reactivity, hence, both models are surface active in their state of maximal 
performance, since for large reactivities the interior becomes rather inactive due to lacking 
particles of type A. The efficiency R is plotted in Fig. [12] a) for L = 2. Similar to the case 
of total selectivity (Fig. El a)) R undergoes a maximum for some intermediate c, however, 
we observe a shift to smaller reactivities. Plotting the maximal efficiency for increasing iV 
(Fig. EDb)) proves the fact that the MTC effect is present for any L. 
Ack 
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VI. CONCLUSION 



Our simulations and analytical results describe the MTC effect quantitatively over a wide 
range of parameters for a model of alternating a- and /3-channels and novel channel topology. 
The output current scales linearly with the number of channels and, hence, the efficiency ratio 
(compared with a topologically and structurally similar reference system without MTC) does 
not sharply decrease with increasing system size. This new insight is encouraging as it allows 
for the MTC effect in experimental settings not only for nanoscale systems. This suggests 
that MTC may enhance significantly the effective reactivity in zeolitic particles with suitable 
binary channel systems and thus may be of practical relevance in applications. Introducing 
a model of size selectivity where reactant as well as product particles are allowed to enter 
a-channels reveals only minor differences in the observed behaviour. For the important case 
of small channel segments, size selectivity seems to be slightly more efficient compared to 
total selectivity. 
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